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TECHNICAL NOTE 3598 

METHOD AND TABLES FOR DETERMINING THE TIME RESPONSE TO A 
UNIT IMPULSE FROM FREQUENCY-RESPONSE DATA AND FOR 
DETERMINING THE FOURIER TRANSFORM OF A 
FUNCTION OF TIME 

By Carl R. Huss and James J-; Donegan 
SUMMARY 


A simple and rapid method is presented for the determination of the 
time response to a unit impulse from frequency-response data and for 
evaluating the Fournier tr ans form of a function of time. Both methods are 
applicable to linear functions for which Fourier, transforms exist. Tables 
are presented which greatly reduce the time required to perform the compu- 
tations. Procedures for performing the calculations by use of the tables 
are illustrated and outlined step by step. 


INTRODUCTION 


The direct and inverse Fourier transforms are at present being used 
extensively in the analysis of flight test data. The direct transform 
provides the means of determining the frequency response of a system from 
a transient response to an input. The inverse transform provides the 
means of determining the time response to either a unit impulse or an 
arbitrary input from the frequency response of a system. In both cases 
the transfer of data between the time plane and the frequency plane is 
accomplished without knowledge of the transfer function relating the 
input and output. 

Some recent experience in analyzing and extending flight measurements 
obtained with a swept-wlng bomber has led to the establishment of a simpli- 
fied method similar to that of reference 1 of determining the time response 
to a unit impulse from frequency-response data and of determining the 
Fourier transform of a function of time. The simplification of the methods 
results from the use of prepared tables. Since these methods and tables 
should be of help to others engaged in the analysis of linear systems, they 
are included in the present paper along with explanatory material of the 
mathematical principles involved, examples showing the application of the 
method and tables, and several comparisons illustrating the accuracy. of 
the procedures. 
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SYMBOLS 


modulus of Fourier transform 


F(icu) 


Fourier transform. 


f(t)e ^dt 


H(irn) 


function of time 


frequency response 


time response to a unit impulse (total unless with, subscript n) 

imaginary component 

interval 

real component 

amplitude — 


t or r 


time, sec 


output 


normalizing parameter, ^pt or ^kx> 


incremental values 


phase angle of Fourier transform, deg 
frequency, radians /sec 


input 


Subscripts : 


refers to particular rectangle 

output 

input 


TnaYiTnirm 


Absolute value of a term is denoted by 


max 


NACA TN 3598 


3 


METHOD 


The methods presented in this paper are applicable to linear systems 
for vhich Fourier transforms exist. The methods are presented in the 
following order: (l) determination of the time response to a unit impulse 

from frequency-response data and (2) determination, of the Fourier trans- 
form of a function of time. Also included as appendix A is a procedure 
for numerically evaluating Duhamel 's integral with an integrating matrix 
based on the Newton-Cotes quadrature formulas given in reference 2 which 
permit the determination of the response of a system to an arbitrary input 
when the response to a unit impulse is known. 


Determination of Time Response to a Unit Impulse 
From Frequency-Response Data 


As indicated by George F. Floyd in reference 1, the time response to 
a unit impulse h(t) can be related to the frequency response by 



cos cut do* 


( 1 ) 


where R [il(ia))] is the real part of the frequency response. 

The integration of equation (l) is sometimes difficult since an 
analytic expression for R [n(ia>)] is often not known and, when one is 
known, it can be very complicated . Numerical values of R[k(ia>)], how- 
ever, are easily obtained from the amplitude ratio and cosine of the 
phase angle of the frequency response and, therefore, the problem lends 
itself well to some method of numerical integration. 

The numerical method employed by Floyd to solve equation (l) was to 
approximate the R[H(ia>)J with straight lines. This type of approximation 
reduces equation (l) for each straight-line segment to the simple form of 
sin x/x for which tables existed. The total response is then the sum of 
the individual response for each straight-line segment. This method gives 
good results for simple systems ; however, when applied to a system with 
several high-frequency modes, it requires considerable time. 

The method suggested in this paper is to approximate the R[H(ia))J 
with a staircase type of function having equal frequency intervals and of 
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such height that the area under each step of the staircase function 
equals the area under that portion of the R [H(im)J curve within the 
interval. 

Equation (l) can then he written as 

^(■t) = fill h n(t) (2) 

n=l 

where 

P«% ' 

h^t) = I RjH(ito)|cos out cko (3) 

and n, o^, and are shown in figure 1(a) and, for n=l, = 0. 

In equation (3 ) , h^t) represents the contribution to the response to 

the unit impulse of the area of the interval between m Q _ 1 and of 

\ 

the staircase function. By substituting the constant amplitude r n for 
rJh(14o)J in equation ( 3 ) and performing the integration, the individual 
contribution becomes 

hn(t) = r n f sin^pt^cos j(2n - l)^ptj (4) 

where r n and £o) are defined as in figure 1(a). 

Examples of the h n (t) for steps of unit amplitude and a frequency 

interval of 1 radian per second of various locations on the frequency 
scale are shown in figure 2. 

Equation (4) can be normalized by letting 


z = 

2 


(5) 
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Then, 


h n (t) = AMT, 


sin z cos( 2 n - l)z 
z 


( 6 ) 


By substituting equation ( 6 ) into equation (2), the total response is 
given by 


CO . 

■u/ 4 .\ 2 * \ si 11 z cos( 2 n - l)z 

h(t) = > r n 3 >- 


n=l 


( 7 ) 


where the time values corresponding to the z -values for which the calcu- 
lations are made depend on the chosen. 

By tabulating the function s f n . z . c . 9 s /^ . n ~ L ) Z . for various values 

z 

of z and n, the response to the unit impulse can be easily and rapidly 
obtained for a chosen frequency interval tea. Such a tabulation is pre- 
sented as. table I for values of z from 0 to 2-5 and values of n from 
1 to 50 . 

The following procedure is suggested in computing the time response to 
a unit impulse from frequency-response data. The steps of the procedure 
axe outlihed and Illustrated by the use of tables I and II and figure 3- 

(1) Choose the frequency interval to be used. The choice of ZVn 
depends on such factors as the accuracy desired, the shape, and frequency 

range of the R [H(ia))J . In general, best results are obtained if the 
interval is chosen so that the staircase function adequately represents 
the R jH(ico)]. (That is, within the intervals, the R[H(ia>)] is as con- 
stant as practicable). For the example shown in figure 3> a of 1 

was chosen. 

(2) Fit the staircase function to R[H(ia))] by using the chosen Zku 
as shown in figure 3> la this case 20 intervals were used. In order to 
obtain accurate results, rJh( 1 o))J must either be zero or approach zero 
within the number of intervals used. 

(3) Read and record the amplitude of each interval as shown in 
column (g) of table II. 
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(4) Read h^t) for each interval at the desired values of z from 
table I and record as shown in column (3) of table II. The example is 

solved for z = 0.5 which, for dtu = 1 radian per second, corresponds to 
time of 1.0 second. 

( 5 ) Multiply the amplitude and h n (t) together for each individual 
interval as shown in col umn (4) of table II. 

(6) Sum the products of the amplitude and h n for the individual 
intervals at the desired values of z. This total is shown as the summa- 
tion of column (?) in table II. 

( 7 ) Multiply the summation by 2 to obtain the time response 

to a unit impulse. This product is shown at the bottom of table II and 
the result is shown in figure ^ at t = 1.0 second. 

Examples of applications of the method are presented in figures 3 
and 4 along with comparisons of the results with the exact time responses . 
In figure 3 "the method has been applied to a heavily damped linear system 
with one low-frequency mode defined by the transfer function 

^ = . In figure 4 the method has been applied to the system 

5 s 2 + 6s + 10 

of figure 3 with the addition of two higher frequency modes. This system 
was defined by the transfer function 



1 100 225 

s 2 + 6s + 10 s 2 + 0.4s + 100 s 2 + 0.2s + 225 


Determination of the Fourier Transform 
of a Function of Time 

As in the determination of the time response to an impulse from 
frequency-response data, the determination of the Fourier transform of a 
time function involves the fitting of a staircase function, in this case 
with constant time intervals, to the time function so that the area of 
the individual time intervals equals the area under the represented por- 
tion of the function. The real part and the imaginary part of the Fourier 
transform for' the individual intervals are then dete rmin ed and the results 
for the various intervals are s umme d to give the real part and imaginary 
part of the total Fourier transform of a time function. The modulus and 
phase angle of the Fourier transform are then easily obtained. 
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The Fourier transform of a time function is 


Too 

F(ico) = f(t) e" icut dt 

Jo 


which for a staircase function can he written 


where 


F(iiu) 


y~ F n (im) 

n=l 


F n (i«) 



r n e’"' i ' m ^dt 


( 8 ) 


(9) 


where r Q _, t nJ and t DL _ 1 are shown in figure 1(h) and for n = 1, ■ 

tn_i = 0. 

The expression F n (io>) represents the contribution of the individual 
intervals to the total Fourier transform. 

The real part of equation (9) is 


R n [F(ico)J = r n | sin^cojcos j^(2n - 1 )^cdJ 
and the imaginary part is 

I n [F(ico)] = -r n | sin^ujsinj^(2n - l)^cu 


( 10 ) 


( 11 ) 


where r n and At are shown in figure 1(h). 
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The modulus of the Fourier transform for an individual interval is 


An(“> =J{*nf * (^f- 


and the phase angle is 


4> n (m) = -tan" 1 

K n 


( 12 ) 


(15) 


Examples of the Fourier transform for a time interval of 0.1 second 
at various locations in the time plane are shown vectorially in figure 5 • 

The real part and imaginary part for the total Fourier transform of 
a time function is simply the sum of the real part and imaginary part of 
the individual intervals . 


R[F(iw)] =^R n [F(im)j (1*0 

and 

ijV(ito)] = ^ I n [F(io>)] ( 15 ) 


In a manner similar to the time response to a unit impulse the real and 
imaginary part (eqs. (10) and (ll)) can he normalized hy letting 


z 



( 16 ) 


Then the real part (eq.. (l0)) becomes 


R n [F(ia))] = Atr n g - ln - z - c ° g . ^ - Q - - -^ ) z - 


(17) 



V 
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and tlie imaginary part (eq.. (11)) becomes 


In[?(M] - ( 18 ) 


The total real part is then 


’ R (F(iai)] 



( 19 ) 


and the total imaginary pant is 


l[p(icD)l = -At \ r ^ z sin(£n - l)z 
• z 


( 20 ) 


The total modiolus of the Fourier transform is 


A(to) = 


y^R [F(im)]| 2 + |l[F(ia>)]p 


( 21 ) 


and the total phase angle is 


$(co) = -tan -1 ' (22) 

R[F(icD)] 


Again, as in the case of the time response to a unit impulse, the 
calculation of the Fourier transform of a time function can be simplified 

by the use of tabulated values of the function s - n — cos (2°- -J-) z _ Rr)f ^ 

z 

s_ in z s i n( 2 n l) _ z ^ or var j_ 0 us values of n and z. It should be noted 
z 
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that the function. sin z cos ^ 2n ~ . - L ) z i s the same function that was 

z 

tabulated to calculate the time response to a unit impulse . The real 
part of the Fourier tr ans form can, therefore, he calculated, from the 
same table. It is only necessary, then, that the additional func- 
tion a - j ~ rl ' — he tabulated. This tabulation has been done 
z 

for values of z from 0 to 1.0 and values of n from 1 to 50 and is 
presented as table III. 

In order to interpolate between the z-values of tables I and III, 
plots can be made for each interval by using the values given in the 
table . 

Fourier tr ans forms of lengthy time functions using a small time 
increment At may be computed by using the tables in conjunction with 
the shifting theorem. This theorem states that, if the Fourier trans- 
form F(icn) of a time function f(t) exists, the Fourier transform of 

f(t - t) is e~ iair F(iu>) . Therefore, a time function extending beyond 
the limits of the table in time may be partitioned into segments , the 
first extending from 0 to t^, the second extending from to t 2 , 

the third extending from r 2 to t^, and so forth. The total Fourier 

transform of the time . function is simply the vectorial sum of the Fourier 
transforms of each segment. The Fourier transform of the first segment 
is obtained directly from the tables. The Fourier transform of each 
succeeding segment is obtained by assuming the origin to be at r and 
then multiplying the results obtained from application of the tables by 

e- icwr . 


The following procedure is suggested for computing the Fourier trans- 
form of a function of time. The steps of the procedure are outlined and 
illustrated by the use of table I, table III, table IV, and figure 6. 

(l) Choose the time interval to be used. The At used depends on 
the highest frequency desired from the Fourier transform. If it is 
assumed that one complete oscillation in the time plane can be defined by 
no less than 4 points,, the mayinrum frequency for which accurate results 
can be expected will be 


^ YnftY 


2jt 

3 At 


2.094 

At 


( 23 ) 
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When the desired mpnri mirm frequency does not dictate the interval to. he 
used, then the interval should he chosen so that the time function is as 
constant as practicable within the interval. For the example shown in 
figure 6, a At of 0.1 second was chosen. 

(2) Fit staircase functions to the time function as shown ih figure 6. 
In order to obtain accurate results at the low frequencies, the function 
of time must damp to practically zero within the number of intervals used. 
For this example 25 Intervals were used, 

(3) Read and record amplitude of each interval as shown in column (2) 
of table IV. 

(4-) Read the values for the calculation of the real part of the 
Fourier transform for each interval at the desired z-values from table I 

and record in column ( 3 ) of table IV. The example is solved for z = 0.25 
which for At = 0.1 corresponds to to = 5 radians per second. 

(5) Read the values for the calculation of the imaginary part of the 
Fourier transform for each interval at the desired z-values from table III 
and record as shown in column (5) of table IV. The example is solved for 
z = 0.25 which for At = 0.1 corresponds to co = 5 radians per second. 

( 6 ) Multiply the amplitude and the values for the real part and the 
amplitude and the values for the imaginary part together for each indi- 
vidual interval as shown in columns ( 5 ) and ( 6 ) , respectively, of table IV. 

( 7 ) Sum the products of the amplitude and the values for the real part 
and the amplitude and the values for the imaginary part. These totals are 

shown in table IV as the summation of columns ( 5 ) and ( 6 ), respectively. 

(8) Use equation ( 19 ) to compute the real part of the total Fourier 
transform and equation (20) to compute the imaginary part of the total 
Fourier transform. These computations are demonstrated at the bottom of 
table IV. 

(9) Use equation (21) to obtain the amplitude of the Fourier transform 
and equation (22) to obtain the phase angle of the Fourier transform. 

This calculation is demonstrated at the bottom of table IV and the results 
obtained are shown in figure 6 at o = 5 radians per second. 

An example of the application of this method to a time function and 
a comparison of the results with the exact amplitude and phase angle is 
shown in figure 6. The time function shown in figure 6 is the same time 
response to a unit impulse as is shown in figure 3- In figure 6 the 
Fourier transform of this time function was obtained and in figure 3 the 
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real part of the Fourier transform was used to obtain the time function 
of the inverse Fourier transform. 


DISCUSSION 


As was mentioned in the "Method" section, in order for the methods 
of this paper to give correct results, the staircase function must ade- 
quately represent the function being fitted. This will be the case if 
(l) the interval used is not too large, ( 2 ) the function being fitted 
approaches zero within the number of intervals used, and (3)- the area 
under the staircase function within the interval equals the area under 
the represented portion of the function being fitted. 

When the time response to a unit impulse is obtained, too large an 
interval may mask some of the modes present in the frequency response; 
thus, in the time response to the unit impulse, such modes would not be 
properly represented either as to frequency or amplitude. If the real 
part of the frequency response does not approach zero within the number 
of intervals used, large errors in h(t) will be obtained near zero 
time. This result is demons t rate d in figure 3 where the error at zero 
time could be reduced to 0.0080 by using 50 intervals and fitting the 
curve to a frequency of 50 radians per second. If the area of the stair- 
case function does not equal the area of the real part of the frequency 
response for each interval, then, of course, inaccuracies will be obtained 
throughout the calculations. 

An idea of the accuracy of the computations involved In determining 
the time response to a unit impulse by the method and tables of this 
paper can be obtained from figure 3 for a simple system and from figure 4 
for a system with multiple modes. The overall accuracy would, of course, 
be increased either by vising a smaller frequency interval or by vising 
more intervals . 

The total time required to obtain response to a unit impulse such 
as that shown in figures 3 and 4 by application of the method presented 
herein and vising a desk calculator is about 1.5 hours. 

In the case of obtaining the Fourier transforms of a function of 
time, the interval chosen determines the maxi mum frequency that can be 
expected to be determ in ed accurately from the method. The larger the 
interval, the lower will be the maximum frequency obtained. If the time 
function does not approach zero within the number of intervals used, an 
end correction must be added in order to obtain the correct Fourier 
transform. The greatest inaccuracies, however, will be in the very low 
frequency range near zero radians per second. As in the previous case 
if the area under the staircase function does not equal the area of the 
time function for each interval, inaccuracies are obtained throughout the 
calculations . 
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The accuracy of the calculations involved in evaluating the Fourier 
transform by the method and tables of this paper can be seen in figure 6. 
As before, the overall accuracy could, of course, be increased by using 
a smaller time interval but, in view of the comparison shown, it would 
not seem to be necessary for the function shown. 


The time required if a desk calculator and the method and tables of 
this paper are used to evaluate the Fourier transform of a time function 
such as shown in figure 6 is about 2.0 hours. 


It should be noted that in some cases the values of z given in 
the tables may not correspond to sufficient time or frequency values to 
define accurately the h(t) or F(ido). The values of time and frequency 
corresponding to the z-values in the tables, as are shown by equations ( 5 ) 

and ( 16 ), depend on the interval chosen. The functions - ^i 1 — - ° - ~1 ~U Z - 


and 


sin z sin(2n - l)z 


can be determined for values of z not given in 


the tables by interpolation (that is, by plotting the numbers given in the 
table for each interval and reading the value of the function of the 
desired z from the plots ) . 


The time response to an impulse when used with Duhamel's integral 
permits the time response to any arbitrary input to be computed; hence, 
the response of a given system to inputs beyond those normally or safely 
experienced by the system may be obtained. This procedure is outlined in 
appendix A. As an example of this procedure, the time response to a tri- 
angular input was calculated by using Duhamel's integral with the response 
to a unit impulse obtained by the method of this paper for the multiple -mode 
system shown in figure 4. The comparison of the calculated time response 
with the exact time response for the arbitrary input shown in figure 7 
serves to demonstrate the procedure and the accuracy. The time required 
to go through the operations with a desk-type computing machine is about 
1 hour. 


The evaluation of the Fourier transform of a time function is also 
useful for many purposes. One of the most important uses is the determina- 
tion of the frequency response of input -output data without knowledge of 
the transfer function involved. The frequency response is simply the 
Fourier transform of the output divided by the Fourier transform of the 
input. The comparison shown in figure 8(c) of a frequency response 
obtained by using the method of this paper to evaluate the Fourier trans- 
form of the input and output time histories shown in figures 8(a) and 8(b) 
with the exact frequency response indicates that, for frequencies below the 
bottoming frequency of the input (about 12 radians per second), the results 
are very good. As would be expected the results above the bottoming fre- 
quency were poor since, at frequencies greater than the bottoming fre- 
quency of the input, the frequency content of the input is very small and 
in such cases the values obtained for a frequency response are unreliable. 
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The time required to calculate the frequency response shown in figure 8(c) 
with a desk-type computing machine is approximately 2.5 hours. 


CONCLUDING REMABKS 


A simple and rapid method is presented for the determination of the 
time response to a unit impulse from frequency-response data. Also 
included is a similar method of evaluating the Fourier transform of a 
function of time . The methods are applicable to linear functions for 
which Fourier transforms exist and give accurate results within the limi- 
tations of the values chosen for the variables. 

Tables are presented which facilitate the necessary computations of 
the methods . These tables reduce the amount of labor and thus the time 
required for the necessary computations . Equations axe presented which 
can be used to obtain greater accuracy or to extend the range of the 
variables. Suggested procedures for performing the necessary computations 
of the methods are illustrated and outlined step by step. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 22, 1955* 
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APPENDIX 

NUMERICAL EVALUATION OF DUHAMEL INTEGRAL 


Duhamel’s integral may be stated as 


x(t) = 6(t) A(O) + F h(t) 5(t - t) dT (Al) 

Jo 


where A(0) is the indicial admittance at t = 0, 5(t) is the time 

history of the input, and h(t) is the response of the system to a unit 
impulse. For the usual aircraft inputs, A(0) =0 and 


x(t) 



h(t) 5(t - t) dT 


(A2) 


This equation may he expressed numerically as 


x ( t i) = At |Z C ij 6 ( t i-j) j = 0, 1, 2, 3, . . . i (A3) 


where 

At incremental time interval 

8(t) time history of the arbitrary input 

h(t) time response to a unit impulse input 

Cj_j integrating matrix given in table V 

Equation (A3) represents i + 1 equations. In order to illustrate 
its use, consider an arbitrary input S(t.) and a system with a known 
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response to a unit impulse h(t) . The response of this system to 
the arbitary input 8(t) is then computed for a time interval 
At = 0.05 second as follows: 


x(0) = 0 


x(0.05) = 0.05 jo.5h(0) 8(0.05) + 0.5h(0.05) 8(0)] 


x(O.l) = 0.05 [o.333331i(0) 8(0.1) + 

1-33333311(0.05) 8(0.05) + 0.3333311(0.1) 6(0)] 


x(0.15) = 0.05[0.3T5h(0) 8(0.15) + 1.125h(0.05) 8(0.1) + 
1.125h(0.1) 8(0.05) + 0.375h(0.15) 6(0)3 

x(0.2) = 0.05 |o.311h(0) 6(0.2) + 

1.422h(0.05 ) 8(0.15) + 0.5331i(0.l) 8(0.1) + 
1.422h(0.15) 8(0.05) '+ 0.3111i(0.2) 8(0)] 


and so forth. The integrating matrix as presented herein has been 
to be very accurate. 


(A4) 


(A5) 


(a6) 


(A7) 

found 
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Value of function 9 10 T — ccs ^ ~ 1 ) z a t values of z of - 


1.0000 0.9999 0.9997 

1.0000 .9995 -9981 

1.0000 .9988 .9949 

1.0000 .9975 -9901 

1.0000 .9959 -9937 

1.0000 .9939 -9758 

1.0000 .9916 .9663 

1.0000 .9888 .9552 

1.0000 .9856 .9427 

10I 1.0000 .9820 .9286 



31 1.0000 

32 1.0000 

33 1.0000 

34 1.0000 

35 1.0000 

36 1.0000 

37 1.0000 

38 1.0000 

39 1.0000 

40 1.0000 

41 1.0000 

42 1.0000 

43 1.0000 

44 1.0000 

45 1.0000 

46 1.0000 

47 1.0000 

48 1.0000 

49 1.0000 

50 1.0000 


0.04 

0.05 

0.06 

0.07 

0.08 

0.09 

0.9989 

0.9984 

0.9976 

0.9967 

0.9957 

O.9947 

• 9935 

.9884 

.9832 

■ 9772 

•9702 

•9625 

.9798 

.9687 

.9547 

-9386 

.9201 

■ 8992 

.9608 

•9390 

.9426 

.8816 

.8464 

.8069 

.9356 

.9001 

.8572 

..8074 

.7510 

.6886 

.9044 

.8522 

■7895 

•7173 

.6365 

• 5480 

.8675 

• 7958 

.7105 

.6132 

• 5056 

•3897 

.8251 

•7314 

.6212 

.4972 

-3619 

.2187 

.777V 

.6597 

.5231 

-3713 

.2090 

.0407 

•7246 

.5814 

-4173 

.2383 

.0507 

-.1386 

.6673 

.4974 

.3056 

.1005 

— 1089 

-3134 

.6056 

.4083 

.1895 

-.0391 

-2657 

-.4781 

•5401 

-3152 

.0707 

-.1781 

—415 6 

—6274 

.4712 

.2189 

-.0492 

-3135 

-5551 

-7563 

•3992 

.1205 

-.1683 

-4428 

-.6803 

-.8609 

-3247 

.0208 

—2850 

-5634 

—7881 

-.9376 

.2480 

-0791 

-3977 

-.6732 

-8758 

-.9840 

.1699 

-.1782 

-5045 

-7696 

-.9412 

-9987 

.0907 

-2755 

-.6041 

-.8510 

-.9825 

-.9862 

.0107 

-3700 

-6952 

-9158 

-9987 

-9317 

-.0691 

-.4609 

-7761 

-9625 

-9893 

-8523 

-.1486 

-.5471 

-.8459 

-9906 1 

-9547 

-7452 

-.2271 

-.6279 

-9036 

-9992 

-.8958 

-.6142 

-.3042 

-.7024 

-.9481 

-.9882 

-.8139 

-4635 

-3794 

-7899 

-9791 

-9579! 

-7112 

-2974 

-.4520 

-.8292 

-9961 

-90891 

-5904 

-.1219 

-5218 

-.8817 

-9986 

-.8421 1 

-4545 

-0575 

-5883 

-9253 

-9869 

-.7588 

-.3070 

-2351 

-.6510 

-9563 

—9608 

-.6607 

-.1516 

.4051 

-.7096 

-9825 

-.9211 

-5496! 

.0076 

.5620 

-.7636 

-9946 

-8681 

-.4278 

.1666 

.7007 

-.8127 

-9978 

—8026 

-2976 

-3214 

.8167 

-.85 66 

-9950 

-7255 

-.1616 

.4680 

•9064 

-.8950 

-9794 

—6380 

-.0224 

.6026 

.9668 

-.9278 

-9502 

--5413 

.1172 

•7218 

•9960 

-9546 

-9174 

-.4368 

-2546 

.8226 

-9930 

-9752 

—8740 

-.3261 

-3869 

• 9024 

.9580 

-9897 

-.8196 

—2107 

.5117 

• 9591 

.8918 

-9978 

-.7616 

-.0922 

.626k | 

•9913j 

■7970 

-9995 

-.6884 

.0276 

.7289 

• 9982 

.6763 

-9929 

-.6156 

.1470 

.8171 

•9796 

-5339 

-.9804 

-5314 

.2642 

.8894 

-9360 

•5741 

-9665 

-.4483 

.3778 

-9442 

.8684 

.2023 

-9430 

-3539 

.4858 

.9806 

-7787 

.0240 

-9134 

-2568 

.5868 

-9978 

.6691 

-1552 

-.8780 

-.1616 

.6795 

-995 4 

• 5424 

-3293 

-8370 

-.0623 

-7623 

.9736 

.4019 

—4928 

-.7908 

.0376 

.8342 

-9328 

.2510 

-.6404 

-7393 

.1371 

.8941 

.8736 

-0938 

-.7672 

-.6832 

.2353 

.9411 

-7974 

-.0658 

-8693 
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TABLE I.- TABULATION QF THE FUNCTION B ^. 1 coa(2° - l? z ._ Contlnued 

z 






Value of funct 

. _ sin z cos(2n - l)z 

values 

of z 

Of - 






z 



O.15 

0.16 

0.17 

0.18 

0.19 

0.20 

0.21 

0.22 

0.23 

0.24 

0.25 

0.26 

0.27 

0.28 

0.29 

1 

0.9851 

0.9830 

0.9808 

0.9785 

0.9761 

0.9736 

0.9709 

0.9680 

0.9651 

0.9620 

0.9588 

0.9557 

0.9521 

0.9486 

0.9448 

a 

.0971 

.8832 

.8685 

.8531 

.8368 

.8198 

.8021 

-7836 

.7644 

•7446 

•7241 

.7030 

.6812 

.6588 

.6358 

3 

.7290 

.6937 

.6567 

.6182 

.5782 

.5367 

.4940 

.4499 

.4049 

.3588 

.3120 

.2645 

.2164 

.1678 

.1188 

if 

.4957 

.4338 

.3899 

■ 3041 

.2371 

.1688 

• 0999 

.0306 

-0389 

-.1080 

-.1764 

-2439 

-3100 

-3746 

-4370 

5 

.2182 

.1298 

.0406 

-.0489 

-.1380 

-.2257 

-.3115 

-3947 

-.4745 

-.5504 

-.6216 

-6879 

-.7481 

—8023 

—8499 

6 

-.0788 

-1873 

-.2934 

-.3958 

-4932 

-5846 

-.6688 

’-7446 

-.8115 

-.8684 

-9147 

-9499 

-.9734 

-9851 

-9848 

7 

-.3688 

-4854 

-5937 

-.6918 

-.7781 

-.8512 

-.9098 

-9528 

-9798 

-.9902 

-9838 

-.9609 

-.9216 

-8669 

-7977 

8 

-.6258 

-7342 

-.8260 

-8992 

-9521 

-9834 

-9927 

-9795 

-9444 

-.8882 

-.8120 

-7178 

-.6076 

-.4839 

-3495 

9 

-.8270 

-.9085 

-.9639 

-99 13 

-9901 

-.9604 

-.9030 

-8195 

-7127 

-5854 

-.4415 

-.2850 

-.1206 

.0470 

.2129 

10 

-.9542 

-9905 

-9913 

-.9563 

-.8869 

-.7857 

-.6564 

-5035 

-3327 

-1503 

.0372 

.2231 

.4007 

.5635 

.7057 

11 

-.9962 

-9720 

-9052 

-7987 

-6572 

-.4870 

-.2956 

-.0916 

.1163 

.3187 

.5068 

.6724 

.8079 

.9078 

.9677 

12 

-9493 

-.8548 

-7155 

-.5387 

-3337 

-.1114 

.1164 

.3378 

-5412 

.7157 

.8522 

•9437 

.9852 

.9749 

.9131 

13 

-8175 

-.6508 

-.4440 

-2097 

.0374 

.2818 

• 5084 

.7030 

.8536 

.9510 

.9891 

.9658 

.8822 

•7441 

• 5599 

Ilf 

-.6127 

-.3808 

-.1215 

.1463 

.4032 

.6305 

.8118 

■ 9340 

.9885 

.9713 

.8837 

•7324 

.5281 

.2860 

.0257 

15 

-3532 

-.0720 

.2149 

.4835 

.7114 

.8796 

•9742 

.9872 

.9180 

-7721 

.5620 

•3055 

•0237 

-2595 

-.5204 

18 

-.0621 

.2440 

.5266 

.7587 

.9182 

.98 99 

.9673 

.8523 

.6566 

.3984 

.1027 

-.2022 

-.4874 

-7256 

-8943 

17 

.2345 

• 5353 

.7780 

.9366 

.9939 

-9439 

•7922 

• 5552 

.2586 

-0653 

-3817 

-.6565 

-8599 

-9702 

-9756 

18 

• 5102 

.7723 

.9405 

.9944 

.9279 

-7489 

.4794 

.1522 

-.1930 

-5143 

-.7727 

-.9372 

-.9876 

-.9183 

-7379 

19 

.7403 

.9308 

.9952 

.9248 

-7295 

.4356 

JO833 

-2798 

—6045 

-.8470 

-9745 

-9701 

-.8343 

-5859 

-.2588 

20 

.9042 

.9948 

.9360 

-7366 

.4270 

.0536 

-3273 

-6585 

-8907 

-9883 

-.9577 

-.7466 

-.4436 

-0745 

-3049 

21 

.9875 

.9578 

.7696 

-4539 

.0636 

-3369 

-.6810 

-.9118 

-9912 

-.9062 

-.6714 

-3256 

.0734 

.4596 

.7689 

22 

.9824 

.8235 

.5150 

.1131 

-.3089 

-.6742 

—9164 

-9913 

-.8858 

-.6194 

-.2406 

.1814 

.5695 

.8534 

-9315 

23 

.8897 

.6058 

.2016 

-.2422 

-.6374 

-9051 

-9924 

-.8820 

-5964 

-1925 

.2491 

.6404 

.9035 

.9864 

.8729 

2 k 

.7175 

.3264 

-1349 

-.5665 

-.8748 

-9931 

-8959 

-.6047 

-.1829 

■2779 

.6777 

•9302 

.9804 

.8181 

.4789 

25 

.4811 

.0139 

-.4560 

-.8182 

-9874 

-.9242 

-.6438 

-.2122 

.2686 

.6855 

-9405 

.9741 

-7785 

•3999 

-0717 

26 

.2005 

-2999 

-7249 

-9650 

-9592 

-.7096 

-2796 

.2207 

.6642 

•9381 

.9723 

.7604 

■3547 

-.1405 

-5989 

27 

-.0961 

-5834 

-.9108 

-.9880 

-7942 

-.3820 

.1330 

.6116 

.921B 

.9787 

.7667 

.3457 

-.1698 

-.6380 

-9302 

2o 

-3843 

-.8076 

-9924 

-.8844 

-.5158 

.0044 

.5226 

.8860 

.9877 

•7982 

. -3759 

-.1602 

-.6461 

-9405 

-9573 

29 

-.6378 

-.9498 

-.9604 

-.6674 

-.1638 

• 3911 

.8214 

.9916 

.8484 

.4372 

-.1119 

-.6240 

-9384 

-9558 

-6713 

30 

-.8358 

-.9956 

-.8185 

-.3648 

.2115 

.7151 

•9773 

.9083 

.5326 

-.0226 

-.5685 

—9226 

-9637 

-.6791 

-.1656 

31 

-.9586 

-.9402 

-.5828 

-0155 

.5566 

.9256 

.9634 

.6520 

.1061 

-4773 

-8897 

-9775 

-.7147 

-1949 

•3941 

32 

-9978 

-.7895 

-.2803 

.3358 

.8224 

•9946 

.7820 

.2715 

-.3426 

-.8241 

-9883 

--7739 

-.2624 

• 3488 

.8250 

33 

-.9422 

-5586 

.0541 

.6441 

.9708 

.8992 

.4648 

-.1607 

-.7198 

-9847 

-.8463 

-3657 

.2647 

-7859 

.9860 

34 

-.8078 

-2709 

.3825 

.8698 

.9807 

.6706 

.0 666 

-5623 

-.9475 

-.9227 

-.5021 

.1391 

.7163 

.9831 

.8245 

35 

-.6006 

.0442 

.6669 

.98if0 

.8508 

■ 3292 

-3430 

-.8568 

-9782 

-.6522 

-.0251 

.6073 

•9642 

.8798 

• 5934 

36 

-3358 

.3549 

.8751 

.9720 

.5993 

-.0624 

-6930 

-.9880 

-8055 

1.2342 

.4483 

■ 9148 

• .9376 

.5078 

-.1663 

37 

-.0464 

.6295 

.9831 

.8355 

.2624 

-4435 

-.9226 

-9311 

-.4655 

.2366 

.8167 

.9805 

.6442 

-.0193 

-.6718 

38 

.2524 

.8398 

-9784 

.5917 

-.1119 

-.7561 

-9918 

-6967 

-.0285 

.6”>40 

-9847 

.7870 

.1675 

-5406 

-.9574 

39 

-5239 

.9656 

.8618 

.2722 

-4703 

-94 69 

-.8887 

-5297 

.4143 

•9235 

•9127 

.3855 

-.3568 

-.8967 

-.9300 

40 

.7509 

■ 99 30 

.6466 

-.0825 

-.7615 

-.9891 

-.6310 

• .1002 

•7710 

.9844 

.6111 

-.1180 

-•7797 

-.9788 

-5983 

Ifl 

.9106 

.9195 

■3573 

-.4262 

-.9442 

-8749 

-2637 

.5110 

.9674 

.8227 

.1686 

-5903 

-.9807 

-.7620 

-0709 

h 2 

.9908 

•7527 

.0271 

-7154 

-.9921 

-.6227 

■ 1495 

.8245 

.9627 

.4838 

-.3210 

-.9064 

-.9025 

-.3124 

.4796 

43 

-9788 

' .5095 

-.3062 

-.9129 

-.8985 

-.2764 


.9808 

■7579 

.0202 

-7299 

-9831 

-.5675 

■2327 

-8733 

44 

.8790 

.2146 

-.6045 

-.9934 

-.6767 

.1224 


.9503 

•3955 

-.4393 

-9593 

-7997 

-.0710 

.7067 

.9814 

45 

.7063 

-.1022 

-.8336 

-.9466 

.-3583 

.4951 

.9801 

.7389 

—0492 

-7995 

-9589 

-.4051 

.4456 

-9647 

.7684 

46 

.4665 

-.4085 

-.9671 

-7784 

.0111 

• 7910 

• 9592 

.3866 

-4835 

-9791 

-.7181 

.0967 

.8556 

.9281 

•3042 

47 

.1854 

-.6735 

-.9900 

-.5103 

•3790 

.9625 

■7716 

-0393 

-.8174 

-9373 

-.3038 

.5730 

.9876 

.6079 

-.2596 

48 

-.1122 

-.8699 

-.8996 

-.1768 

.6928 

.9821 

.4499 

-4577 

-.9813 

-.6838 

.1850 

•8977 

-8587 

.1021 

-7384 

49 

-3997 

-.9782 

-.7061 

.1792 

.9078 

.8466 


-.7890 

-.9412 

-2756 

.6284 

.9851 

.4854 

-4350 

-9757 

50 

-.6516 

-.9870 

-.4318 

.5124 

•9932 

•5775 

Bess 

-.9699 

-.7054 

.1948 

.9180 

.8121 

—0261 

-8391 

-.8939 
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NACA TN 3598 



Value of function sin z cog(2n - l) z a ^. of z of - 


1.05 1.10 1.15 1-20 1.25 1.50 1-35 1-40 1-45 1.50 1-55 I 1.60 | 1.65 | 1.70 | 1-75 


.4231 -5742 .6834 .7456 .7589 .7239 .6455 .5307 .3888 

.3987 .12 44 -.1546 -.4034 —5929 —7024 -.7226 —6549 — 5124 

-.8258 -.7205 -.4773 -.1507 .1911 .4800 .6611 .7035 -6061 



-3324 -.4163 
.6742 .0406 
•7484 .3467 
•5244 -.6345 
.0913 -7412 

-3771 -8354 
.6954 .3478 
.7380 .0394 
.1939 -4153 
.2503 .6723 

.4190 -.7370 

.7154 .5906 

.7270 -2752 
.4490 -.1190 
.0076 .4791 

.5602 —7021 
.7298 -7241 
-7104 -5389 
.4065 .1994 
.0603 .1972 


.7207 -.3031 
.6750 .4984 
-5002 —6358 
.2285 .7005 
.0865 -.6836 

■3849 .5878 
.6095 -.4242 
.7171 -2115 
.6871 .0256 
-5253 -.2598 

.2627 .4639 
.0502 -.6144 
.3536 .6940 
.5891 -6933 
.7116 .6125 

.6975 -.4610 
.5496 .2561 
.2965 -.0217 
.0139 -.2179 
•3214 .4273 


.3893 .6627 
.2439 -.6639 
.0830 .6518 
.0819 -.6258 
.2422 .5891 

.3884 -.5398 
-5U.9 -4795 
.6059 -.4096 
.6645 .3316 
.6846 -.24 69 


1.55 

1.60 

1.65 

0.0134 

-0.0182 

-0.0478 

-.0402 

.0347 

.1422 

.0670 

-.0908 

-2330 

1-0938 

.1268 

.3181 

.1200 

-.1624 

-.3951 

-.1463 

.1972 

.4623 

.1723 

-.2315 

-5179 

-.1980 

.2649 

.5605 

.2234 

-2975 

-5890 

-.2483 

.3291 

.6029 

.2728 

-.3600 

-.6013 

-.2964 

.3887 

.5853 

-3205 

-.4166 

-.5343 

-.3435 

.4430 

•5095 

.3659 

-.4680 

-.4518 

-.3876 

.4913 

.3829 

.4088 

-.5130 

-.3044 

-.4292 

•5330 

.2182 

.4488 

--5510 

-.1266 

-4677 

-5673 

.0319 

.4857 

-.5816 

.0637 

-.5030 

• 5939 

-1579 

.5194 

-.6042 

.2477 

-5348 

.6125 

-3315 

•5493 

-.6186 

.4070 

-5629 

.6226 

-.4724 

.5755 

-6245 

•5259 

-5871 

.6243 

-.5662 

•5977 

-.6220 

• 5924 

-.6073 

.6175 

-6057 

.6158 

-.6109 

.6000 

-6233 

.6023 

-.58H 

.6296 

-5915 

.5478 

-.6348 

.3788 

-.5007 

.6391 

-.5641 

.4411 

-.6422 

.5475 

-3704 

.6441 

-.5291 

.2905 

-.6449 

-5087 

-2033 

.6447 

-.4867 

.1110 

-6433 

.4630 

-.0160 

.6408 

-.4377 

-0795 

-.6372 

.4110 

.1730 

.6325 

-.3828 

-.2622 

-.6267 

• 3533 

.3448 

.6198 

-3227 

-.4187 

-.6118 

■ 2909 

.4822 

.6028 

-.2581 

-5336 

-5928 

.2245 

.5716 

.5817 

-.1901 

-5953 

-.5696 

.1551 

.6041 
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Vklua of function 


Bin z cos(2n - l)z 


at values of z of - 


1.80 

1.85 

I.90 

1.95 

2.0 

2.05 

2.10 

2.15 

2.20 

2.25 

2.30 

2-35 

2.40 

2.45 

-0.1229 

-0.1432 

-0.1610 

-0.1764 

-0.1892 

-0.1996 

-0.2075 

-0.2151 

-0.2163 

-0.2172 

-0.2160 

-0.2128 

-0.2075 

-0.2005 

.3434 

-3861 

.4158 

.4324 

.4385 

.4290 

.4111 

.3839 

.3492 

.3088 

.2644 

.2180 

.1712 

.1257 

-.4929 


-.4987 

-.4514 

-3815 

-2937 

-1955 

-.0948 

.0016 

.0870 

.1567 

.2074 

.2389 

.2474 

-5407 

-.4768 

.4818 

.3700 

.2231 

.0622 

-.0914 

-2193 

-.3080 

-3502 

-5455 

-.2996 

-.2232 

-.1296 

-.0334 

-.3056 

-.0886 

.1276 

.3002 

.3988 

.4106 

.5418 

.2137 

.0586 

-.0895 

-.2018 

-.2602 

-2599 

•3143 

.0365 

-.2298 

-.4083 

-.4548 

-3671 

-.1832 

.0342 

.2189 

• 3208 

.3197 

.2282 

.0841 

-.0635 

-0873 

.2438 

.4322 

.4852 

.2908 

.0232 

-.2309 

-.3890 

-.3482 

-■1939 

.0178 

.1962 

.2749 

.2382 

-.1580 

-.4498 

-4855 

-.2671 

.0701 

.3404 

.4097 

.2816 

-.0052 

-2393 

-■3257 

-2329 

-0358 

.1516 

.3706 

.5193 

.5158 

-.0779 

-.3860 

-.4144 

-.1704 

• 1595 

■ 3512 

.2947 

.0549 

-.1904 

—2812 

-1796 

-5070 

-4307 

-.0141 

•3795 

.4342 

.1382 

-.2422 

-3893 

-.2110 

.1148 

• 3114 

.2378 

-.0132 

-.2186 

.5382 

.2119 

-2934 

-4735 

-.1819 

.2580 

.4083 

.1528 

-.2215 

-3431 

-.1247 

.1845 

.2789 

.0981 

-.4586 

.0717 

.4784 

.3081 

-.1964 

-.4327 

-.1581 

.2868 

.3471 

.0298 

-.2834 

-.2424 

.0820 

.2552 

.2842 

-3335 

-.4833 

.0263 

.4382 

.2396 

-2532 

-.3668 

.0081 

.3305 

.1883 

-.1785 

-.2680 

-.0029 

-.0512 

.4939 

.2545 

-.3462 

-.3770 

.1574 

.4064 

.0271 

-3521 

—1692 

.2411 

.2468 

-.1088 

-2563 

-.1924 

-.5044 

.0806 

.4764 

.0542 

-.4204 

-.1453 

.3450 

.2033 

-2591 

-.2424 

.1724 

.2490 

-0927 

.3963 

.3815 

-3505 

-3454 

.3061 

.3280 

-.2640 

-.3036 

.2240 

.2784 

-.1888 

-.2511 

.1524 

.2217 

-.5183 

-.1089 

.4938 

.0232 

-4545 

.0456 

-4o4l 

-.1016 

-.3461 

.1417 

.2843 

-.1661 

-.2223 

.1754 

.9334 

-.1768 

-4507 

.3089 

.2879 

-3785 

-1323 

.3850 

-.0113 

-3382 

.1230 

.2552 

-1913 

-1583 

-.4383 

.4088 

.1875 

-.4737 

.0780 

.3895 

-.2744 

-.2071 

-3530 

.0008 

-.3118 

.1598 

.1888 

-2336 

.2527 

-.5168 

.1340 

.3788 

-.3899 

-.0693 

.4013 

—2191 

-2057 

-3578 

-.0531 

-.2592 

.2244 

.0891 

-.0149 

•4875 

-3993 

-.0763 

.4317 

-.3098 

-.1191 

.3826 

-.2266 

-.1433 

• 5237 

-1534 

-.1496 

-2594 

-2259 

-.2764 

.4980 

-.2681 

-1745 

-4255 

-.2845 

-0877 

.3450 

-2774 

-.0196 

.2630 

-.2505 

.0276 

.4201 

.0012 

-3883 

.4654 

-.2036 

-.1794 

■ 3981 

-3123 

.0148 

.2602 

-3194 

.1469 

.1057 

—2491 

-.5276 

.2743 

.1163 

-.4078 

.4407 

-2193 

-1059 

■ 3381 

-3539 

.1677 

.0912 

-.2866 

.2690 

-.1206 

.5261 

-.4664 

.2044 

.1265 

-3725 

.4315 

-2943 

.0413 

.2030 

-3310 

.2989 

-.1403 

-0587 

.2042 

-.4160 

.3169 

-.'4396 

.2241 

.0463 

-.2768 

• 3945 

-3712 

.2291 

-.0282 

-1583 

.2701 

-2793 

.1967 

.2200 

-.4103 

.4911 

-.4518 

.3120 

-U.33 

-0925 

.2582 

-3438 

-3429 

-.2634 

.1338 

.0098 

-.1308 

.0214 

.1792 

-3372 

.4319 

-.4541 

.4070 

-.5038 

.1658 

-.0178 

-.1164 

.2175 

-2734 

.2810 

-2455 

-.2584 

.1065 

.0424 

-1753 

.2817 

-3548 

• 3903 

-.3891 

.3548 

-2938 

.2147 

-.1268 

-0394 

.0392 

.4421 

-3598 

.2701 

-1774 

.0859 

.0007 

-.0789 

.1481 

-.2003 

.2402 

-2655 

•2765 

—2741 

.2801 

-•5545 

.9038 

-.4897 

.4329 

-.3940 

-5538 

-3129 

.2720 

-.2316 

.1925 

-1551 

-1199 

-0873 

.0578 

-5164 

-.4947 

-4729 

-.4511 

.4292 

-4075 

.3858 

-.3642 

• 5427 

-.3214 

.3003 

-2795 

.2588 

-2386 

-3918 


-.2784 

.2220 

-.1671 

.1146 

-0653 

.0199 

.0210 

-.0570 

.0877 

-.1130 

.1326 

-.1488 

.1863 

-.0741 

-0325 

.1287 

-.2107 

■2757 

-3217 

•3482 

-3558 

.3454 

-.3200 

.2823 

-.2356 

.1838 

.0577 

—20 98 

-3298 

-.4089 

.4426 

-.4318 

.3808 


.1976 

-.0886 

-.0160 

.1060 

-1738 

.2154 

-.2898 

.4297 

-.4892 

.4850 

-.3679 

.2205 

-.0517 

-.1085 

.2342 

— 5080 

.3031 

—2849 

.2053 

-1055 

.4620 

-5192 

.4442 

—2662 

.0583 

.1781 

-3301 

.3860 

-.3416 

.2185 

—0566 

-.0989 

.2098 

-2539 

-5389 

.4510 

-.2134 

-.0785 

.3178 

-.4252 

•3754 

-.2010 

-.0242 

• 2159 

-.3108 

.2874 

-.1685 

.0087 

•5045 

-.2458 

-.1065 

.3802 

-4537 

.3108 

-0379 

-.2249 

• 3564 

-.3096 

.1264 

.0918 

-2392 

.2572 

-3659 

-.0340 

.3820 

-4734 

.2754 

■0679 

-3382 

.3812 

-.1948 

-.0854 

.2825 

-.2897 

.12 68 

.0872 

.1518 

.3036 

-.4978 

•3072 

•0937 

-.3889 

• 3895 

-.0807 

-.2366 

.3456 

-.1898 

-.0846 

.2614 

—2247 

■0936 

-.4809 

.4054 

.0274 

-3979 

.3792 

-.0242 

-.3166 

• 5405 

-.0603 

-.2399 

.2918 

-.0869 

-.1710 

-3198 

.3121 

-.1436 

-3470 

.4285 

-.0470 

-5445 

.3344 

.0275 

-.3201 

.2438 

-0774 

-2756 

.1609 

.4799 

-3877 

-1783 

.4764 

-1598 

-.3234 

• 3832 

.0485 

-•3572 

.1952 

• 1853 

-.2937 

.0327 

.2310 

-.5409 

.1456 

.4258 

-3447 

—2178 

.4208 

-.0104 

-3733 

.1921 

.2378 

-.2851 

-.0701 

.2813 

-.0747 

.4902 

.l4o8 

-4950 

.0240 

.4443 

-.1587 

-3531 

.2508 

.2391 

-2955 

-.1213 

.2954 

.0165 

-2589 

-.3383 

-.3844 

.3574 

.3098 

-3831 

—2384 

.3386 

.1723 

-.3391 

-.1132 

• 3124 

.0628 

—2784 

-.0218 

.1166 

.51 12 

-0735 

-4738 

.0303 

.4321 

•0055 

-388 9 

-0307 

• 3435 

.0513 

-.2970 

-0653 

•2507 

.1292 

-4839 

-.2459 

.3781 

• 5234 

-.2591 

-.3800 

• 1395 

• 3579 

-0315 

-.3238 

-0554 

.2670 

.1154 

-.3483 

.3076 

• 4595 

-.0751 

-.4532 

-.1348 

.5495 

.2771 

-.1893 

-.5300 

.0214 

.2983 

.1120 

-2077 
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TABLE II.- ILLUSTRATION OF SUGGESTED PROCEDURE TO BE 

USED TO CALCULATE THE TIME RESPONSE TO A UNIT 

IMPULSE FROM FREQUENCY -RESPONSE DATA 

(Values are from fig. 3 ancL table I at z = O .5 which 
for a = 1 corresponds to a time of 1.0 sec] 



I 


0 

n 

r n 

knOt) 

0 

X 

0 

1 

0.0937 

0.8415 

0.0788 

2 

. 05^9 

.0678 

.0037 

3 

.0157 

-.7682 

-.0121 

4 

-.0050 

-.8979 

.0045 

5 

-.0123 

-.2021 

.0025 

6 

-.0135 

.6795 

-.0092 

7 

-.0126 

,9364 

-.0118 

8 

-.0111 

• 3324 

-.0037 

9 

-.0096 

--5773 

.0055 

10 

-.OO83 

-.9562 • 

.0079 

11 

-.0071 

-.4560 

.0032 

12 

-.0061 

.4634 

-.0028 

13 

-.0053 

.9568 

-.0051 

14 

-,oo46 

.5705 

-.0026 

15 

-.0042 

-.3403 

.0014 

16 

-.0037 

-.9382 

.0035 

17 

-.0033 

-.6735 

.0022 

18 

-.0029 

.2104 

-.0006 

19 

-.0027 

.9007 

-.0024 

20 

-.0024 

.7631 

-.0018 

• 



y~ 0 = 0.0611 


h(l.O) = — Zxjd 



0 = |xlx (0.06ll) 


0.0390 
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TABLE m.- TABULATION CF TEE FUNCTION sl B - g Bln J 2P ~ 1 ^ Z USED TO CALCULATE 
THE IMAGINARY PART OF THE FOURIER TRANSFORM OF A FUNCTION OF TIME 




Value of the function a ^ n 

z aln(2n - l)z 
z 

at values of z 

of - 



0 

0.05 

0.10 

O.15 

0.20 

0.25 

0.30 

0-35 

0.40 

0.45 

0.50 

1 

0 

0.0500 

0.0997 

0.1489 

0.1974 

0.2448 

0.2911 

0-3359 

0.3791 

0.4204 

0.4597 

2 

0 

.1494 

.2950 

•4333 

.5609 

.6746 

•7716 

.8498 

.9074 

.9431 

.9565 

3 

0 

-2473 

.4786 

.6791 

.8359 

-9391 

.9826- 

.96k) 

.8853 

.7521 

-5739 

4 

0 

■ 3428 

.6431 

.8642 

.9789 

•9738 

.8503 

.6248 

.3261 

-.0081 

-3364 

5 

0 

.4348 

.7820 

-9721 

.9674 

.7700 

.4210 

-.0082 

-.4308 

—7622 

-9373 

6 

0 

.5225 

.8897 

.9932 

.8031 

■ 3777 

-1554 

-.6374 

-.9264 

-9394 

-.6765 

7 

0 

.6049 

.9619 

-9255 

-5121 

-.1071 

-6775 

-9668 

-.8601 

-.4058 

.2063 

8 

0 

.6814 

•9958 

.7752 

.1402 

-5656 

-9629 

-.8415 

—2720 

.4350 

.8994 

9 

0 

.7510 

.9900 

.5558 

-2538 

-.8857 

-.9120 

-.3204 

.4810 

.9466 

.7657 

10 

0 

.8131 

.9447 

.2864 

-.6078 

-.9889 

-5425 

•3514 

•9423 

.7418 

-.0720 

11 

0 

.8671 

.8617 

-.0084 

-8658 

-.8500 

.0166 

.8579 

.8320 

-.0244 

-8435 

12 

0 

•9124 

.7444 

-.3024 

-9871 

-5030 

.5698 

.9609 

.2170 

-7721 

-.8395 

13 

0 

.9486 

•5975 

-5694 

-9526 

-.0329 

•9240 

.6121 

-5296 

-9355 

-0637 

lk 

0 

•9753 

.4267 

-7856 

-7676 

.4454 

•9554 

-.0247 

-9550 

-.3910 

.7708 

13 

0 

■ 9923 

.2388 

-.9316 

-.4615 

.8145 

.6531 

-6498 

-.8010 

.4494 

.8965 

16 

0 

■ 99 94 

.0415 

-9943 

-.0825 

•9843 

.1226 

-9693 

-.1612 

.9497 

.1980 

17 

0 

.9965 

-1575 

-9683 

-3095 

.9130 

-4507 

-.8330 

.5764 

.7313 

-.6825 

18 

0 

.9856 

-.3502 

-8557 

.6526 

.6182 

-8665 

-.3048 

-9643 

-.0406 

-9355 

19 

0 

.9609 

-.5289 

-6668 

.8927 

.1721 

-9797 

■ 3667 

-7674 

-7817 

-.3284 

20 

0 

.9286 

-.6866 

-.4182 

-9919 

-.3162 

-.7506 

.8657 

.1049 

-9313 

.5806 

21 

0 

.8870 

-.8169 

-1323 

•9345 

-7270 

-2593 

.9576 

-.6212 

-3761 

.9558 

22 

0 

.8366 

-.9146 

.1654 

•7295 

-9599 

.3226 

-5991 

-9705 

.4638 

.4523 

23 

0 

.7778 

-.9759 

.4484 

.4094 

-9578 

.7918 

-.0411 

-7311 

.9526 

-.4671 

2k 

0 

.7112 

-.9982 

.6913 

.0246 

--7211 

.9844 

-6621 

-.0483 

.7205 

-9571 

25 

0 

-6375 

-.9808 

.8724 

-.3640 

-.3079 

•8331 

-9716 

.6639 

-.0568 

-.5671 

26 

0 

•5575 

-.9242 

-9757 

-.6952 

.1807 

•3909 

-.8242 

-9733 

-7912 

.3442 

27 

0 

.4718 

-.8309 

.9917 

-.91 66 

.6250 

-.1880 

-.2891 

.6924 

-.9268 

-9391 

20 

0 

.3815 

-.7044 

.9192 

-9933 

.9164 

-.7012 

.3819 

-.0086 

-.3610 

.6706 

29 

0 

.2874 

-.5498 

-7646 

-9131 

-9833 

-9694 

.8733 

-.7042 

-4779 

-.2145 

30 

0 

.1903 

-3733 

■5417 

-6889 

.8095 

-.8990 

• 9539 

-9728 

• 9552 

-9023 

31 

0 

.0915 

-.1819 

.2704 

-3559 

.4376 

-5145 

.5861 

-.6513 

.7096 

-.7606 

32 

0 

-.0084 

.0168 

-.0251 

•0334 

— 0416 

.0496 

-0576 

-0653 

-.0730 

.0805 

33 

0 

-.1082 

.2147 

-3183 

•4173 

-.5105 

.5965 

-.6741 

.7423 

-.8003 

-8475 

3k 

0 

-.2068 

.4041 

-5831 

•7354 

-.8544 

-9350 

-•9735 

-9690 

—9220 

.8349 

35 

0 

-3034 

•5774 

-7957 

•9374 

-.9892 

.9469 

-.8152 

.6079 

-3459 

.0543 

36 

0 

-3969 

.7278 

-.9374 

.9913 

-.8817 

• 6279 

-2733 

-1220 

.4920 

-7763 

37 

0 

-.4866 

.8490 

-.9955 

.8888 

-.5584 

.0897 

.3969 

—7778 

.9576 

-8934 

38 

0 

-.5714 

.9384 

-.9642 

• 6459 

-.0984 

-4798 

.8807 

-.9618 

.6978 

-.1898 

39 

0 

-.6503 

.9865 

-.8471 

• 3012 

.3854 

-8818 

•9501 

-5625 

-.0902 

.6884 

40 

0 

-.7229 

•9972 

-6543 

-.0913 

•7755 

-9757 

•5727 

.1782 

-.8099 

•9337 

41 

0 

-.7882 

.9683 

-.4030 

-.4692 

•9753 

-.7288 

-.0740 

.8106 

-.9170 

•3205 

42 

0 

-.8457 

.9007 

-1157 

-7731 

•9364 

-2273 

-.6859 

■ 9512 

-.3307 

-.5873 

43 

0 

-.8946 

• 7971 

.1819 

-.9550 

.6682 

•3537 

-•9753 

.5143 

.5059 

-.9552 

44 

0 

-9346 

.6619 

•4633 

-9859 

.2364 

.8110 

-.8059 

-.2346 

■ 9596 

-.4448 

45 

0 

-.9654 

■5001 

•7032 

-8614 

-2532 

.9851 

-2576 

-.8412 

.6872 

.4745 

46 

0 

-.9864 

.3186 

.8803 

-.6008 

-.6809 

.8151 

.4120 

-9378 

-.1054 

.9576 

47 

0 

-•9977 

.1243 

.9789 

-.2453 

-.9418 

.3603 

.8877 

-.4660 

-8l8l 

.5603 

48 

0 

-9989 

-.0750 

■9099 

.1488 

-.9722 

-.2204 

•9457 

.2884 

-.9118 

-.3521 

49 

0 

-.9901 

-.2712 

.9126 

.5195 

-.7645 

-.7240 

•5585 

.8680 

-3154 

-.9408 

50 

0 

-.9715 

-.4567 

•7537 

.8081 

-.3696 

-9748 

-.0914 

•9209 

.5196 

-.66^5 
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TABLE IH.- TABULATION OF THE FUNCTION 


sin 8 S ln(2n_^l)z _ Concluded 
z 



Value of 

the function Eln z - 1 ) 2 

z 

at values of 

z of - 



0.55 

0.60 

0.65 

0.70 

0-75 

0.80 

O.85 

0-90 

0.95 

1.00 

1 

0.4967 

0.5314 

0.5655 

0.5929 

0.6195 

0.6433 

0.6640 

0.6818 

0.6965 

O.708I 

2 

.9474 

.9165 

.8650 

.7944 

.7071 

.6057 

.4929 

-3720 

.2462 

.1188 

3 

.3627 

.1328 

-.1007 

-.3228 

-5195 

-.6786 

-.7911 

-.8508 

-.8556 

-.8069 

• k 

-.6183 

-.8202 

-.9188 

-.9042 

-.7806 

-.5661 

-.2891 

.0146 

.3071 

-5528 

5 

-9237 

'•7273 

-.3909 

.0155 

.4090 

-7117 

.8656 

.8442 

.6571 

.3468 

6 

-.21 96 

.2932 

-7097 

.9094 

.8385 

• 5245 

.0661 

-3982 

-.7319 

-.8415 

7 

•7244 

-9397 

.7706 

.2937 

-.2904 

-.7423 

-.8826 

-.6632 

-.1838 

•3535 

8 

.8768 

.3878 

-2975 

-.8096 

-.8796 

-.4812 

.1614 

.6996 

.8508 

-5472 

9 

.0710 

-.6586 

-.9297 

-.5689 

.1659 

•7704 

.8410 

-3453 

-.3662 

-.8090 

10 

-. 8 i 24 

-.8652 

-.1999 

.6162 

.9031 

.4362 

-.3781 

-.8565 

-.6140 

.1261 

11 

-.8080 

.0316 

.8228 

.7784 

-.0382 

--7959 

-.7436 

-0439 

.7632 

•7040 

12 

.0794 

.8881 

.6401 

-.3516 

-.9085 

-3897 

.5897 

.8366 

.1205 

-.7121 

13 

.8800 

.6120 

-.4803 

-.8979 

-.0904 

.8186 

.5968 

-.4240 

-.8411 

-.1114 

14 

.7190 

-.4446 

-8971 

.0464 

.8957 

-5419 

-.7235 

-.6438 

.4234 

.8048 

15 

-.2278 

-.9342 

.0004 

-9137 

.2171 

-.83 86 

-.4104 

.7166 

.5674 

-5584 

16 

-.9256 

-.2325 

-8973 

.2642 

-.8650 

-.2929 

.8292 

.3183 

-.7902 

-3393 

17 

-.612.9 

.7657 

•4797 

-.8239 

-3395 

-8557 

■ 1967 

-.8612 

-.0556 

.8414 

18 

.3704 

.7874 

-.6407 

-.5443 

.8169 

.2430 

-8799 

.0739 

.8265 

-.3610 

19 

.9480 

-.1942 

-.8224 

.6388 

.4551 

-.8699 

.0309 

.8278 

-.4788 

-.5416 

20 

.4896 

-.9288 

.2006 

.7615 

-7526 

-.1922 

.8720 

-.4500 

-.5176 

.8110 

21 

-5030 

-.4780 

-9297 

-.3800 

-.5608 

.8811 

-.2556 

-6239 

.8128 

-.1354 

22 

-9467 

.5823 

.2968 

-.8904 

.6737 

-1399 

-.8065 

.7328 

-.0079 

-.7001 

23 

-3550 

.9001 

-7710 

.0772 

.6561 

-.8894 

.4626 

.2909 

-.8077 

.7160 

24 

.6246 

.0700 

-.7092 

.9170 

-.5809 

-.0888 

.6873 

-.8650 

-5301 

.1041 

25 

-9217 

-.8494 

-3915 

-2336 

-.7388 

.8945 

-.6397 

.1021 

.4649 

-.8026 

26 

.2114 

-6855 

.9186 

-.8375 

.4757 

.0566 

-5225 

.8186 

-.8307 

-5639 

27 

-7297 

-3525 

.0991 

-.5190 

.8062 

-.89 66 

-7743 

-.4741 

.0722 

-3332 

28 

-8735 

.9410 

-.8656 

.6607 

-.3616 

.0158 

-3229 

-.6032 

.7840 

-.8413 

29 

-.0627 

.3285 

-.5628 

.7436 

-8574 

•8957 

-.8576 

.7482 

--5791 

.3669 

30 

.8167 

-7029 

.5641 

-.4078 

.2403 

-.0681 

-.1019 

.2632 

-.4096 

-5359 

31 

.8030 

-.8383 

.8647 

-.8823 

.8914 

-.8918 

.8838 

-.8678 

.8440 

-.8130 

32 

-.0887 

-0947 

-.1015 

.1080 

-.1142 

.1202 

-.1258 

.1311 

-.1561 

.1407 

33 

-.8834 

.9070 

-9190 

.9190 

-.9073 

.8847 

-.8514 

.8082 

-.7559 

.6958 

34 

-7135 

.5626 

-3902 

.2044 

-.0142 

-.1718 

• 3452 

-4983 

.6249 

-7198 

35 

•2359 

-.4993 

-7102 

-.8495 

.9055 

-.8747 

-7825 

-.5818 

.3519 

-.0967 

36 

.9274 

-.9244 

.7701 

-.4932 

-1423 

.2229 

-•5417 

.7626 

-.8523 

.8004 


.6054 

-.1707 

-.2982 

.6819 

-.8854 

.8616 

-.6229 

-2352 

-1992 

-5694 


-3781 

.8007 

-9297 

.7250 

-.2676 

-2732 

•7022 

-.8695 

.7236 

-.3264 

39 

-.9484 

-7510 

-1992 

-.4354 

.8475 

-.8458 

.4420 

.1599 

-.6671 

.8411 

40 

-.4824 

-.2564 

.8231 

-.8730 

.3875 

• 3226 

-.8160 

.7969 

-2923 

-3735 

41 

.5109 

-.9369 

.6396 

.1386 

-.7927 

.8268 

-2317 

—5221 

.8560 

-5301 

42 

.9458 

-.4226 

-.4809 

.9201 

-.4996 

-3708 

.8758 

--5597 

-.2612 

.8148 

43 

.3471 

.6306 

-.8968 

.1742 

-7220 

-.8052 

.0068 

.7763 

-.6872 

-.1480 

44 

-.6309 

.8796 

.0011 

-.8609 

.6018 

-4179 

-8773 

.2069 

.7055 

-.6916 

45 

-9195 

.0068 

.8975 

-.4669 

-.6369 

.7808 

.2201 

-.8703 

.2310 

-7257 

46 

-2033 

-.8747 

.4791 

.7022 

-.6919 

-.4635 

.8206 

.1886 

-.8548 

.0894 

4-7 

• 7351 . 

-.6407 

-.6412 

.7055 

• 5390 

-.7538 

-.4315 

-7847 

-3217 

--7981 

48 

.8702 

.4103 

-.8221 

-.4624 

.7681 

-5075 

-.7094 

-.5451 

.6469 

.5748 

49 

.0544 

.9381 

.2014 

-.8627 

-.4303 

■7242 

.6143 

-5369 

-.7399 

.3196 

50 

-.8209 

.2694 

.9298 

.1691 

-.8290 

-.5498 

• 5511 

.7891 

-.1684 

-.8408 
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TABLE IV.- ILLUSTRATION OF SUCJffiSTED PROCEDURE TO BE USED TO 
CALCULATE THE FOURIER TRANSFORM OF A FUNCTION OF TIME 


^alues are from fig. 6, table I, and table TTT at 
z = 0.25 which for a At = 0.1 corresponds to 
a frequency of 5.0 radians /sec] 


® 

© 

© 

® 

© 

© 

n 

r n 

R[F(iaj)] 

l[F(ie>2 

© x © 

© X ® 

1 

0.0420 

O.9588 

0.2448 

0.0403 

0.0103 

2 

.0950 

■7241 

.6746 

.0688 

.0641 

3 

.H65 

.3120 

• 9391 

.0364 

.1094 

4 

.1200 

-.1764 

.9738 

-.0212 

.1169 

5 

.1120 

-.6216 

.7700 

-.0696 

.0862 

6 

.1000 

-.9147 

• 3777 

-.0915 

.0378 

7 

.0860 

-.9838 

-.1071 

-.0846 

-.0092 

8 

.0720 

-.8120 

-.5656 

-.0585 

-.0407 

9 

.0580 

-.4415 

-.8857 

-.0256 

-.0514 

10 

.0460 

.0372 

-.9889 

.0017 

-.0455 

11 

.0560 

.5068 

-.8500 

.0182 

-.0306 

12 

.0290 

.8522 

-.5030 

.0247 

-.0146 

13 

.0220 

.9891 

-.0329 

.0218 

-.0007 

14 

.0170 

.8837 

.4454 

.0150 

.0076 

15 

.0130 

.5620 

.8145 

.0073 

.0106 

16 

.0100 

.1027 

.9843 

.0010 

.0098 

17 

.OO65 

-.3817 

.9130 

-.0025 

.0059 

18 

.0050 

-.7727 

.6182 

-.0039 

.0031 

19 

.0040 

-9745 

.1721 

-.0039 

.0007 

20 

.0030 

-•9377 

-.3162 

-.0028 

-.0010 

21 

.0020 

-.6714 

-.7270 

-.0013 

-.0015 

22 

.0015 

-.2406 

-9599 

-.0004 

-.0014 

23 

.0010 

.2491 

-.9578 

.0003 

-.0010 

24 

.0005 

.6777 

-.7211 

.0003 

-.0004 

25 

.0002 

-9405 

-3079 

.0002 

-.0001 





Y1 © = -0.1297 

r © = 0.2644 


R[F(i5.0)] - At) (g) = 0.1 X (-0.1297) = -0.01297 (See eq. (19)0 


l[F(i5.03 = “ - 0 - 1 x (0.2644)= -0.02644 (See eq. (20).) 


*(5-0) 


A(5-0) =/R[F(iffi)] 2 + l[F(im)] 2 = 0.02945 (See eq. (21).) 

- -tan -1 = -tan -1 2.0387 or 4> = -U6°08' (See eq. 

r[fM 


( 22 ).) 
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Time, t, sec 

(b) Terms used to obtain the Fourier transform of a function of t ime 

Figure 1.- Illustration of terms used in the development of the 

methods of this paper. 
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Figure 3-~ Comparison of the response to a unit impulse calculated by 
the method of this paper with the exact response for a simple system 

Y "1 

the transfer function of which is £ = . 

® s^ + 6s + 10 


o»|w 


0 .5 1.0 1.5 Z0 2.5 

Time, t, sec 

Figure 4.- Comparison of the response to a unit impulse calculated "by 
the method of this paper with the exact response for a multiple-mode 
system the transfer function of which is 

1 100 22 £ 

s 2 + 6s + 10 s 2 + 0.4s + 100 s 2 + 0.2s + 225 



Time, t, sec 



Figure 5- - Vector representation of the frequency content of an element 
of the staircase function at various locations in the time plane. 
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Frequency, oj, radians/ sec 


Figure 6.— Comparison of the modulus and p has e angle of Hie Fourier 
transform of a time function calculated hy the method of this paper 
with the exact values. 




Response, x Amplitude of input, S 
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Figure 7 .- Comparison of the response to an arbitrary input calculated 
by using the h(t) obtained by the method of this paper vith that 
calculated by using the exact h(t) for the mul tiple-mode system. 
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(a) Modulus and phase angle of input. 

Figure 8 .- Comparison >71111 the exact values of the frequency response 
obtained by application of the method of this paper to an Input 
and output . 


Modulus, A x i$x’ deg Amplitude of output, x 
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